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Derivation of Hydrodynamics for Multi-Relaxation Time Lattice Boltzmann using the 

Moment Approach 
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Department of Physics, North Dakota State University, Fargo, ND 58108, U.S.A. 

A general analysis of the hydro-dynamic limit of multi-relaxation time lattice Boltzmann models is 
presented. We examine multi-relaxation time BGK collision operators that are constructed similarly 
to those for the MRT case, however, without explicitly moving into a moment space representation. 
The corresponding 'moments' are derived as left eigenvectors of said collision operator in velocity 
space. Consequently we can, in a representation independent of the chosen base velocity set, generate 
the conservation equations. We find a significant degree of freedom in the choice of the collision 
matrix and the associated basis which leaves the collision operator invariant. Therefore we can 
^— ^ ■ explain why MRT implementations in the literature reproduce identical hydrodynamics despite 

being based on different orthogonalization relations. 
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o ■ 

0\ ! I. INTRODUCTION 

The lattice Boltzmann (LB) method is continuing to increase in popularity as a simulation method for fluid me- 
chanics for a wide range of applications from turbulence [l[ to complex fluids Q. A key of its success is the simplicity 
of the algorithm. Instead of discretizing the hydrodynamic equations directly the method is based on an underlying 
microscopic model. Historically the method developed from lattice gases [3| where particles move on a lattice and 
collide on lattice points. Because such a lattice gas model locally conserves mass and momentum the macroscopic 
behavior of the system has to be described by the continuitiy and Navier- Stokes equations [3]. The connections be- 
\ tween the microscopic streaming and collision rules and the macroscopic differential equations is established by taking 
the hydrodynamic limit which requires averaging the locally conserved quantities. This reproduces the Boltzmann 
equation [5]. Performing a Taylor expansion on the discrete Boltzmann equation then leads to a PDE representation 
O . of the discrete evolution equation [6fl . 

At this point there are several routes to proceed. Grad_|7| suggests taking moments of the full Boltzmann equation 
, which is a route that has been taken by other groups @. Alternatively one can formally expand the distribution 
function before taking the moments, which is known as the Chapman- Enskog expansion [§[. Either approach will 
lead to identical results to second order: the continuity and Navier-Stokes equations as well as the heat equation for 
thermal systems. The higher order equations are, however, quite different. Here neither approach has been particul arly 
successful as the Navier-Stokes level equations appear to be appropriate to length-scales close to molecular scale [Tof . 
There are few attempts to derive higher order hydrodynamic equations in the LB context. One exception are multi- 
phase fluids where higher order spatial derivatives giving rise to surface tension have to be taken into account fllj ). 

The development of the method took a major leap when it was discovered that it is feasible to use a Boltzmann- 
level microscopic model [lH, which removes microscopic noise. This approach is referred to as the lattice Boltzmann 
method. Qian et ai. [l3| found that the approach is simplified considerably when the collision operator is written as 
a single-time BGK expression which relaxes local particle distributions towards the equilibrium distribution. To this 
date this represents the most popular flavor of lattice Boltzmann algorithms employed. 

Shortly after the introduction of the single-time relaxation collision operator d'Humieres realized that one can 
extend the BGK collision with a multi-relaxation time (MRT) approach [14j. This replaces the scalar relaxation time 
with a matrix, which then allows for decoupled relaxation of different stress terms. Thus it decouples the different 
transport coefficients and they no longer need to take their ideal gas values. Up to now deriving hydrodynamic 
equations for multi-relaxation time lattice Boltzmann methods has been almost exclusively achieved by Chapman- 
Enskog like expansions. These expansions typically depend on the specific model [9] although a very general (but 
not very detailed) approach was presented by Junk et al. [l5j |. Here we show that it is straight forward to derive 
the hydrodynamic equations through a moment approach in a model-independent manner. Our approach clarifies 
the structure of the collision operator by showing that the stress moments have to be left eigenvectors of the collision 
matrix. This feature has been obscured by current implementations of multi-relaxation time LB algorithms, which, 
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either to reduce the degrees of freedom in the choice of the collision matrix [l4| or because of additional properties 
required for analysis of the method fill Il7j. choose that the eigenvectors of the collision matrix be orthogonal with 
respect to a specific scalar product. 

It may be initially surprising that the freedom to choose orthogonality with respect to different scalar products 
exists. The underlying reason, however, is simple: because the conserved moments are left invariant by the collision, 
these moments can be mixed into the non-conserved moments without changing the collision operator. With this 
observation we then reconcile the existing multi-relaxation time lattice Boltzmann models with the analysis in this 
paper. 



II. LATTICE BOLTZMANN 



The lattice Boltzmann equation (LBE) is a representation of the Boltzmann transport equation [5| with three levels 
discretization taken into account: time t , position x and velocity v. First LB-methods utilized a two body collision 
operator derived from lattice gas methods (Higuera et al. [U[). Later Qian and D'Humieres realized that the collision 
operator could be significantly simplified using a BGK approach pj| as = (l/r)(/P - /,) where /P is the local 
equilibrium distribution [l3|, [TiJ . In the BGK approximation [l9[ the collision integral is replaced by a relaxation term 
that moves the current distribution /(x, v, t) function towards the equilibrium distribution / (x, v, t) For a general 
collision operator fli the basic LBE can then be written as 

/i(x + v,,t + i)-/ i (x,t) = n i (/i,-.. j N ), (i) 

where the fi are the density functions associated with a discrete set of N base velocity vectors v i; a; is the lattice 
position and t is the discrete time with interval At = 1. The velocities are chosen such that the Vi are lattice vectors. 
Since collisions conserve certain quantities such as mass and momentum we require 

5>r^ = o, ( 2 ) 

i 

where the ip^ are the vectors describing the velocity moments of the conserved quantities. The index c only emphasizes 
that these vectors are associated with conserved quantities. We will encounter non-conserved vectors ip a later in this 
paper. The first quantity that has to be conserved in the collision is the local density which has a corresponding 
vector of ip^ c = li where 1, is simply 1 for every i. Momentum must also be conserved in each spatial direction. 
In three dimensions the corresponding ip vectors are tp,' = ipf' c = Vi y , and i\)\ ,c = Vi^ z . We denote the locally 
conserved quantities as density p and momentum j. They are defined through the vectors tp a ' c as 

£ ^' C h = P> £ = 3a- (3) 

i i 

Throughout this paper greek indices a, /3, 7 will generally denote the range of spatial dimensions {x,y, z} and be 
treated under the Einstein summation convention. Latin indices k are used in the context of vector components 
of the lattice Boltzmann base velocity set and are summed over explicitly. 

Most LB models are used to simulate isothermal hydrodynamics and these models are the focus of this paper. 
Thermal models require the conservation of the additional moment vf, which we do not treat here. In principle, 
however, it should be easy to extend the presented approach to thermal systems and generate the corresponding heat 
equation. 

To recover the continuity and Navier-Stokes equations this local equilibrum distribution needs to match the 
first four velocity moments of the continuum Maxwell-Boltzmann distribution. This distribution is / (v) = 

(27r9) 3 /' 2 ex P ( ^ 2 IT ) wnere l oca l velocity is defined as u = j/p and 8 is the temperature. For thermal models 
we would need to match velocity moments. The first four moments sufficient to derive isothermal hydrodynamics are 

£/° = P, (4) 

i 

^"»/° = P U a = 3a, (5) 

i 

^ViaVipf® = p9 + pU a U/3, (6) 
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ViaVipV^f® = p8 (u a 5i3 7 + Up5 la + U^Safl) + pUallpUy + Q al 3 7 . (7) 
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The tensor quantity Q a p-y is an arbitrary correction term and vanishes in the continuum case. However, the typical 
choice is Qa/37 = — PUg Upv^y which allows us to use a much smaller velocity set. The tradeoff are small Galilean 
invariance problems [20|. Note that the conserved moments of the local equilibrium distribution /P and the distribution 
fi are identical because the collision does not change them, i.e. ^i^t^fi = Y^i'4 > i' C fi- 

Depending on the base velocity set the conditions Eqs. (j3K]) may not uniquely define the equilibrium distribution. 
From several different general arguments it is usually found that the explicit form 



/°(p, u ) = P w i 



1 1 / s2 1 
(H — u.v, H K (u.Vj) u.u 



(8) 



is the optimal choice for the isothermal equilibrium distribution for an appropriate choice of the Wi weight constants 
[2l| although other forms have been used (22| . 

When deriving the hydrodynamic equations from the continuous Boltzmann equation using the single-relaxation 
time approximation leads to a fixed ratio of the transport coefficients such as shear viscosity, bulk viscosity, and 
thermal conductivity Q . In the discrete case of lattice Boltzmann the same hydrodynamic equations can be derived 
with transport coefficients containing a renormalized relaxation time u) = (t — 1/2). For ideal gases the predicted 
ratios agree quite well with the experimentally measured values The form of the hydrodynamic equations apply 
not only to ideal but also non-ideal gases and even fluids. Lattice Boltzmann applications usually consider examples 
from this more general class of systems. In these more general cases, however, the ratios of transport coefficient are 
no longer fixed, and it would be advantageous to write a more flexible collision term that allows for independently 
variable transport coefficients. This was accomplished by D'Humieres |14| by considering a multi-relaxation time 
BGK collision operator of the form 

Wi,--- Jn) = J> y [/>,u,0)-/ J (x,t)], (9) 

3 

where A is a collision matrix. If we choose Ay = Sij/r we recover the single- relaxation time collision operator. 
Another numerical rationale for implementing multi-relaxation time Lattice Boltzmann methods is the improvement 
in stability, particularly for high Reynolds numbers [23]. There are some requirements on the collision matrix to 
ensure mass and momentum conservation in the collision. In the single-relaxation time approach the conservation 
laws were respected because the conserved moments of the local distribution /$ and the local equilibrium distribution 
/P are identical. For the multi-relaxation time collision term Eq. @ requires 

Etf'TMtf /,) o. (10) 

i 3 

These equations will be satisfied if we demand that the scalar product of a conserved quantity vector with the collision 
matrix is a linear combination of conserved quantity vectors, i.e. J2i i^i^-^-ij — c a ipj' c for an arbitrary c a . Note here 
that the only physically relevant quantity is the collision operator f2, not the collision matrix A. While different choices 
for c a will lead to different collision matrices, they will not change the collision operator fl. A convenient choice that 
coincides with the single-relaxation time case sets the conserved moments lj and Vi a to the left-eigenvectors of our 
collision matrix with some eigenvalue: 

J2^ A ik = — Ik, (11) 

v ji Via 1 (1^) 



v i ak i 
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where we used the relaxation times r to denote the inverse eigenvalues of the collision matrix. This choice also allows 
us to ensure that A is invertible which, while not strictly necessary, simplifies the formalism. Clearly, the values of 
l/r p and l/Tj a are entirely arbitrary, meaning that r p and Tj a may not appear in the hydrodynamic equations. 



III. HYDRODYNAMIC LIMIT BY THE MOMENT METHOD 



In this section we present a new approach to obtain the hydrodynamic equations for the multi-relaxation time 
lattice BGK equation. We generalize the moment approach familiar from single-relaxation time methods [2(| to the 
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more general MKT formalism. For the multi-relaxation time collision operator we expand the left hand side of Eq. ([T]) 
to second order: 



(d t + v ia d a ) fi + \ (d t + v ia d a ) (d t + vipdp) fi + 0(d 3 ) = A u (fj ~ fj) 



(13) 



This allows us to write the fi in terms of the /P and higher order derivatives as long as A 1 exists: 

i 

This is important because we can express the equilibrium distributions /P in terms of p and u in Eq. ([5]) but not 
the local distributions ft. Here we have made the assumption that both, spatial and temporal derivatives, are small 
quantities of the same order of magnitude. As a byproduct we see that the conservation equations by virtue of Eq. (|10[) 
and the -0"' c being left-eigenvectors of Ay require 



]T n° £ (A- 1 ) [(d t + v ia d a ) fi] = ]T l(dt + v ia d a ) fi] = 0(d 2 ) 1 



(15) 



which we will use later. Replacing all occurances of fi in Eq. (flU)) with Eq. (JT5J) up to second order we obtain 

(d t + v ja d a ) /P - (dt + v ja d a ) ]T [(A- 1 )^ - 1,5 J (d t + v l0 d p ) /° + 0(d 3 ) = £ A fi (/P - /,) . (16) 

i i 

Because we know ff as a function of p and j this is an equation expressed entirely in terms of our hydrodynamic 
variables, except for the collision term. So far the only requirement on the collision Matrix A is that it be invertible 
and fulfill Eq. (|10p . The general approach now to obtain a conservation equations is to take the inner product of the 
conserved quantity vectors ip a,c with Eq. (jT5J) . The collision term then vanishes, we retain no dependencies on the 
and, after some algebra, we obtain the conservation equations. 



A. The Continuity Equation 



To obtain the continuity equation we take the inner product of tpj' c — lj with Eq. (|16p from the left hand side, i.e 
we just sum over Eq. (|16[) while making use of mass conservation in Eq. (jlOp . We get 

1 



h (dt + v ja d a ) /j 5 - £ lj (d t + v ja d a ) 



j j 
We can rewrite the second order terms as 



(d t + v if} d ) /P + 0(d 6 ) = 0. (17) 



(d t + v i0 d p )f? = O(d 3 ), 



(18) 



d * E 1 J v ^ E 



(dt + vipdp)/? = 0(d 3 ), 



(19) 



where we used that both lj and lj conserved quantity vectors so that we can apply Eq. p"5)) . We are 

left with 



E( 5 * + %A)/° + o(d 3 ) = o, 



(20) 



which using Eq. (UJ) and Eq. ([5]) becomes the continuity equation 

d t p + d a (pu a )+0(d 3 ) = 0. 



(21) 
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B. The Navier-Stokes Equation 



As the Navier Stokes Equation describes the conservation of momentum we take the first order velocity moment of 
Eq. (fTB|) and obtain 



£ v o<* (dt + vjffdp) /° - ^ v 3<x (dt + Vj-,d 7 ) ^ 



(d t + v t pdp)f° + O(d 3 )=0. (22) 



The collision term vanishes according to Eq. (fTOI) . We can rewrite the first of the second order terms as 



^£%,£[(a-%-^ 



(dt+Vipdp)f? = 0(d 3 ), 



(23) 



which vanishes to third order due to Eq. (|15[) much like Eq. (|19[) . To evaluate the remaining gradient term 



^7 VjaVj-y J]] 



(A- 1 ),,; - 



(d t + V dp)f? 



(24) 



we need to know the stress moments Y]j VjaVj^ (A x ) „ — 



of the collision matrix. From the single-relaxation 



time derivation (20| we know that these terms lead to the stress terms in the Navier-Stokes equation we wish to 
obtain. Because we want to distinguish between bulk and shear stress now we separate these into a trace and a 
traceless velocity moment 



D 



A; 



(25) 



The key requirement is now that the trace and the (D — 1) + l) elements of the traceless part are left eigenvectors 
of the collision matrix A. For the trace part we demand 



(26) 



where tb is the bulk relaxation time and for the traceless part we require 



£ 



D 



D 



(27) 



where the shear stress relaxation time t$ is the eigenvalue. These eigenvalue equations for the second order velocity 
moments are the key property of the collision matrix that allows us to recover the Navier-Stokes equation. Because 
of the freedom to choose different eigenvalues for the trace and the traceless part we can obtain independent bulk and 
shear stresses. 

What follows is essentially the same derivation as in the single-relaxation time case [2p| , except that we now have 
two stress terms with associated relaxation times that need to be treated independently. We use the eigenvalue 
equations ([26)) and (f27|) in Eq. (|24|) to replace A -1 with the appropriate eigenvalues. The different velocity moments 
are substituted by the expressions in Eqs. (j4]-[7]) and we replace tb — \ = ujb and t$ — \ — ws- We get 



<V* £ £ v i 



(A-% 



Vj a Vj~ 



(A-),, - ^ 



,£> + 2 



d t [ pusus — + p08 ai ) + 6 — —8 ai dfj {pup) + —dp (pu s u s up + Qssp) 



D 



D 



d t I pu a u 7 - pusus—— + dp [pO {u a 8p 1 + up8 ia + u^Sap) + pU a UpU^ -f- QaP~/) 



D 

, D $ ai pup + ^jj- (pususup + Qssp) 



(28) 
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To treat the second order terms further we need two identities we obtain by looking at the first order terms of Eq. 
Inserting the moments (j^J), ([5]) and ignoring all second order terms we get 

d t ( P u a ) = -dp ( P 9S a p + P u a ufi) + 0{d 2 ). (29) 

Using the continuity equation (|21[) . we obtain the second identity 

pd t u a = -pupdpu a - dpp9S a p + 0(d 2 ). (30) 

These two identities and the continuity equation Eq. (|21[) now replace the time derivatives in Eq. 

dyUJB I - 0S aj dp (pup) - ^jj- [u s dp (pOSps + pupu s ) + u$ {pupdpug + dpp9Sp S )] 



+ 



D + 2 



D 



+ <9 7 w s j - u^dp {p9S a p + pu a up) - u a (pupdpu~ ( + dpp95 1 p) + dp (pu a upu 7 + Q Q /3 7 ) 

+ -jj- [usdp {pOSps + pupup) + u 5 (pupdpug + dpp96p S )} + dppO (u a Sp 7 + up8 ia + u 7 S a p) 



o{d 3 ) 



+ w s 



- ^j=j-dp (pu s u s up + Qssp) 

^Odapd-fU^ + ^dad^Qjss 

2 I 

9d 1 (d-fUa + d a u y ) + dpd 7 Q a p 7 - —ddapd^u^ - —d a d 1 Q 1 ss 



0(d 3 



If we now combine the first order terms Eq. (|2"9")l with the second order terms Eq. 
moment of the LBE (|22|) we find the Navier-Stokes equation 



p9 (d<yU a + d a u-y) — 
1 



(31) 



of the first order velocity 



dt (pu a ) + dp (pu a up) = - d a p9 + daUJBjjOpd^Uy + dyUJs 



(32) 



daUBjjdyQ^ss + d^uj s (dpQap-f - d^^Qyss^j + 0(d 3 ) 



= - a, 



2 
D 



+ d a pd 1 u 1 + djriKdjUa + <9 Q u 7 ) — -^-9 7 u 7 ^ Q7 ] + 0(d 2 Q) + 0(d 3 ), (33) 



where u = ^p9[jB — h) is the bulk and -q — p9(rs — |) the shear viscosity. 

In summary we recover the continuity and Navier-Stokes equations in a similar form as found from multi-relaxation 
time approaches with independently adjustable bulk and shear viscosities provided that three conditions are fulfilled: 

1. The first four velocity moments of the equilibrium distribution are given by Eqs. (|3][7]). 

2. The moments of the conserved quantity vectors and Vk a are not altered in the collision step. 

3. The collision matrix has the left eigenvectors VkaVkp — Vk-yVk-y-jy- and Wfc 7 Wfc 7 . 

Unfortunately none of the published multi-relaxation time lattice Boltzmann methods [lil . fullfil this last require- 
ment. This is because we have some additional freedom in combining the if>£ vectors with vectors from the conserved 
quantities as we will explain below. It is interesting to note that we have constraints up to the third order velocity 
moments for the equilibrium distribution, but only up to second order moments for the collision matrix. 

We should mention that the derivation presented here does not impose any requirements on the extra degrees of 
freedom that are typically present in a lattice-Boltzmann implementation. A DDQQ simulation with a Q component 
base velocity set in D dimensions only requires K = 1 + D + D(D + l)/2 base vectors to reproduce isothermal 
hydrodynamics: 1 for the density, D for the momentum components, and D(D + l)/2 for the stress tensor. Our 
derivation makes no assumptions about the structure of the remaining Q — K 'ghost' or kinetic modes or the choice of 
their corresponding relaxation times. Often the relaxation times for these ghost degrees of freedom are uniformly set 
to 1. In this case all possible choices for ghost eigenvectors of the collision matrix lead to identical collision matrices. 
The choice of ghost modes can influence the performance of the LB method if one wants to make use of the freedom 
to choose arbitrary relaxation times [24]. The introduction of fluctuations to the LBM requires careful treatment of 
the ghost modes and their relaxation times (l7| , particularly in the context of boundary conditio ns [25l . Furthermore 
Adhikari and Succi suggested a duality between conserved quantity vectors and ghost modes [26j as guideline for 
constructing base velocity sets for multi-relaxation time implementation. 
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C. Limited Freedom of Choice of the Eigenvectors 



When we required Eqs. (|26|) and (|27|) we ignored that there is a remaining freedom of choice for the eigenvectors. 
To understand this, let us first remember that the relaxation times for the conserved moments r p and Tj a are entirely 
arbitrary by construction. Because the conserved moments of the /j and ff are identical the collision term simply 
can not alter the values of the conserved quantities, independent of the value of r p and Tj a . This also implies that the 
effect of adding multiples of a conserved mode eigenvector ip°j ,c to any of the eigenvectors will still result in suitable 

eigenvectors. Consider an alternative collision matrix A with a left eigenvector (-0^ + ipfy : 



E 



(34) 



Here 0™ is an eigenvector of the original matrix A -1 . The n indicates that it corresponds to a non-conserved quantity 
and r n is the associated eigenvalue. In contrast ip? is an eigenvector that corresponds to a conserved quantity, i.e. p 
or v a , with the associated eigenvalue r c . Now, terms that depend on r c have to vanish because its value is entirely 
arbitrary. Therefore we will only retain the t™^™ terms in the hydrodynamic equations. The collision matrices A 
and A will lead to identical hydrodynamic equations. To illustrate this we reinvestigate the bulk stress component in 
the second order terms in the Navier-Stokes derivation in Eq. ([25)) for the alternative collision matrix A. We replace 



Vjsvjs with {vjgVjs 



Klj) - Kl s 



and use the aforementioned new collision matrix A 1 and obtain 



pcdt E E i( v tf v i& + Kl i) - K H 



A" 



2 6ji 



f? 



i 3 



A" 1 ) -;<5 



2" 31 



d t (^pu s u s ^- + p96 a ^j + D ^ 2 06 ay d l3 {pup) 



d 7 u B K [d t p + dp (pup)} ^j- + d 7 oj p K [d t p + dp (pup)} 



D 



(35) 



For readibility we omit the correction terms from Eq. ([7]) here as no additional third order velocity 

moments are generated by the lj term in the new bulk viscosity eigenvector. The lj contributions lead to additional 
terms consisting of derivatives of the continuity equation. Since these contributions vanish to third order Eq. (|21[) 
the resulting Navier-Stokes equation remains unaffected. If we decided to add a first order velocity moment to one of 
the non-conserved eigenvectors we would find a Navier-Stokes equation instead of the continuity equation here which 
again vanishes to third order. We are thus free to add any vectors corresponding to our conserved quantities to the 
non-conserved eigenvectors. This is the degree of freedom that allows us to impose orthogonality on the eigenvectors 
with respect to different inner products. 

To recover the approach of d'Humieres we now need to require all of the left eigenvectors of A^ be orthogonal, 
with respect to the inner product J^j^T^? = 3nmN n where N n is the norm of the vector -0™ which need not 



be normalized. The only non-orthogonal left eigenvectors here are lj and vj 



We remedy this by applying a 



Gram-Schmidt orthogonalization procedure to find the new orthogonalized bulk stress 



Vj 1 Vj 1 



N 



-h 



(36) 



and thus recover d'Humieres' basis. In contrast recovering the Benzi approach requires that the eigenvectors obey 
orthogonality with respect to the Hermite norm: J^j ipJ l ipj'Wj = 6 mn M n . Again only one pair of eigenvectors is not 
orthogonal, lj and Vj-yVj-y. We apply the same orthogonalization procedure, however, with the new norm and thus 
obtain 



£ ljv«;j W 

— li 

N 



(37) 



as the orthogonal bulk stress vector. While d'Humieres' and Benzi's approaches lead to different collision matrices 
it is important to note that a practical implementation of the approaches is entirely identical. This is because the 
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eigenvectors only differ by a multiple of lj, which is the density eigenvector and therefore a conserved quantity 
eigenvector. 

Let us assume that we have two collision matrixes A and A and two corresponding sets of left eigenvectors that 
only differ by a conserved quantity vector ip a and ip a = ip a + ij] c . Vectors with the same index a correspond to the 
same physical quantity and will thus correspond to the same time constant r°. The eigenvalue equations are then 

ifj a A = , and Q + 4> c ) A = — (V> a + ijj c ) . (38) 

We know that conserved quantity vectors ip c are left eigenvectors of both A and A and Eq. (fTU)) requires that 

E ^ E ^ (/? - fs) = E # E k » (fi -fi) = °> ( 39 ) 

i 3 i 3 

independent of the actual choice of basis. Now the collision operators il and f) can be defined as 

Q = E A tJ (/? - , and Q = £ A,, (/° - /,•) . (40) 

3 3 

Operators are defined by their action on a basis. Therefore we let 51 act on an arbitrary vector chosen from its own 
left eigenvector basis. Using Eqs. (HQ]) , ([55]). and we get 

E = E < E A « (# - fi) ( 41 ) 

i t j 

iE '•:(/;' -o) ( 42 ) 

= ^E(^ + ^)(/°-/i) ( 43 ) 

= ew + ^<)E a «(#-£) (44) 

= E<E A «(/°-/i) (45) 
= E^' (46) 

Thus we have proven that as long as two different bases differ only by conserved quantity left eigenvectors, the collision 
operators are, in fact, identical. 



IV. SUMMARY 



We presented a new general formulation for the derivation of hydrodynamics. Based on the framework of generalized 
or multi-relaxation time formalism we performed a direct asymptotic expansion to second order of the lattice Boltz- 
mann equation and derived the continuity and Navier-Stokes equations for the isothermal ideal gas. Our approach is 
general in the sense that we do not require specific knowledge of the base velocity set and equilibrium distribution 
function as long as the velocity moments to third order are identical to those of the continuous case and the collision 
does not affect the conserved quantities. We therefore do not require an explicit multi-relaxation time representation 
but instead describe all physically relevant quantities in terms of left eigenvectors of a collision matrix. These left 
eigenvectors can again be described in terms of velocity moments and thus we maintain a representation independent 
of the chosen base velocity set. The eigenvalues of the collision matrix are chosen to be inverse of the relaxation time 
related to the physical quantity in question. Through the relaxation times associated with bulk and shear stress terms 
we then get direct access to the bulk and shear viscosities. 

The derivation illuminates a degree of freedom in the choice of the left eigenvectors. This is rooted in the fact that 
the collision does not alter conserved quantities. Therefore linear combinations of conserved quantity eigenvectors 
can be added to non-conserved moment left eigenvectors without changing the collision operator and by extension 
the hydrodynamic equations. We identify this degree of freedom as source of the validity of multi-relaxation time 
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implementations based on different inner products such as the standard vector and the Hermite norm. In fact, we show 
that for the simple case of isothermal hydrodynamics the collision operators of any two realizations of multi-relaxation 
time Lattice Boltzmann are identical provided they conserve mass and momentum and the appropriate equilibrium 
distribution is chosen. 
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